> 
oo 
o 

00 



O 



Progress on tilted axis cranking covariant density functional 
theory for nuclear magnetic and antimagnetic rotation 

J. Meng, 1 ' 2 * 3 J. Peng, 4 S. Q. Zhang, 1 and P. W. Zhao 1 

1 State Key Laboratory of Nuclear Physics and Technology, 
School of Physics, Peking University, Beijing 100871, China 
2 School of Physics and Nuclear Energy Engineering, 
Beihang University, Beijing 100191, China 
3 Department of Physics, University of Stellenbosch, Stellenbosch, South Africa 



O 

^Department of Physics, Beijing Normal University, Beijing 100875, China 



0\ ' (Dated: January 10, 2013) 

: Abstract 

Magnetic rotation and antimagnetic rotation are exotic rotational phenomena observed in weakly 
deformed or near-spherical nuclei, which are respectively interpreted in terms of the shears mecha- 



nism and two shearslike mechanism. Since their observations, magnetic rotation and antimagnetic 
rotation phenomena have been mainly investigated in the framework of tilted axis cranking based 
on the pairing plus quadrupole model. For the last decades, the covariant density functional theory 
and its extension have been proved to be successful in describing series of nuclear ground-states and 
excited states properties, including the binding energies, radii, single-particle spectra, resonance 
states, halo phenomena, magnetic moments, magnetic rotation, low-lying excitations, shape phase 
transitions, collective rotation and vibrations, etc. This review will mainly focus on the tilted 
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axis cranking covariant density functional theory and its application for the magnetic rotation and 
antimagnetic rotation phenomena. 

PACS numbers: 21.60.Jz, 21.10.-k, 21.10.Rc, 23.20.-g 
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I. INTRODUCTION 



A. Basic concept 

The quantal rotation has been recognized and investigated with a distinguished history for 
one hundred years. In nuclear physics, the rotation features of the nucleus were firstly noted 
by Teller and Wheeler in 1938 (l|. A more complete explanation of the nuclear rotation 
was due to Bohr and Mottelson by pointing out that the rotation was a consequence of 
deformation Q and by developing a model that combined the individual and collective 
motions of the nucleons 

Since then it was generally accepted that the high angular momenta of atomic nuclei 
are connected with the collective rotations with a stable nuclear deformation. The nuclear 
rotational bands are built on the states with substantial quadrupole deformations and char- 
acterized by strong electric quadrupole (E2) transitions, which can be interpreted as the 
coherent collective rotation of many nucleons around an axis perpendicular to the symme- 
try axis of the nuclear density distribution jsj. The study along this line has been at the 
forefront of nuclear structure physics for several decades, as evidenced by the exciting dis- 
coveries such as the baekbending g angular momentum alignment fl Q, superdeformed 
rotational bands [7J, etc. 

In earlier 1990s, however, the rotational-like sequences of strongly enhanced magnetic 
dipole (Ml) transitions were surprisingly observed in several light-mass Pb isotopes, which 
are known to be spherical or near-spherical. Sequentially, this new type of rotational bands 
which have strong Ml and very weak E2 transitions has been well discovered experimentally 
in a number of nearly spherical nuclei (for reviews see Refs. Besides the outstanding- 

magnetic dipole transitions, the intriguing feature here is that the orientation of the rotor is 
not specified by the deformation of the overall density but rather by the current distribution 
induced by specific nucleons moving in high-j orbitals. 



The explanation of such bands was given in terms of the shears mechanism In this 
interpretation, the magnetic dipole vector, which arises from proton particles (holes) and 
neutron holes (particles) in high-j orbitals, rotates around the total angular momentum 
vector. See Fig. [1] for a schematic illustration. At the bandhead, the proton and neutron 
angular momenta are almost perpendicular to each other. Along the bands, energy and an- 
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gular momentum are increased by an alignment of the proton and neutron angular momenta 
along the total angular momentum. Consequently, the orientation of the total angular mo- 
mentum in the intrinsic frame does not change so much and regular rotational bands are 
formed in spite of the fact that the density distribution of the nucleus is almost spherical or 
only weakly deformed. In order to distinguish this kind of rotation from the usual collective 
rotation in well-deformed nuclei (called electric rotation), the name "magnetic rotation" 
(MR) was introduced in Ref. Q, which alludes to the fact that the magnetic moment is the 
order parameter inducing a violation of rotational symmetry and thus causing rotational-like 
structures in the spectrum [3]. This forms an analogy to a ferromagnet where the total 
magnetic moment, the sum of the atomic dipole moments, is the order parameter. 



FIG. 1: A schematic illustration for the spin-coupling scheme of magnetic rotation. For a near- 
spherical nuclei, the coupling of the proton-hole j~ l and neutron-particle j u , each in high-j orbital, 
gives the total angular momentum J. As a result, a large transverse component of the magnetic 
dipole moment vector, n±, rotates around the total angular momentum vector, and creates the 
enhanced Ml transitions. 



1) a AI = 1 sequence of strong magnetic dipole transitions, corresponding to a reduced 
transition probability B(M1) ~ a few [a 2 n ; 

2) weak or absent quadrupole transitions, corresponding to a deformation parameter \f3\ < 
0.15, which combined with strong Ml transitions results in large B(M1)/B(E2) ratios, 




The experimental indicators for magnetic rotation can be summarized as 
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3) a smooth variation in the 7 transition energy with angular momentum; 

4) a substantial dynamic moment of inertia, corresponding to the large ratio of the 
jW/B(E2) > lOOMeV-^efe)" 2 , compared with the values in well-deformed [~10 
MeV _1 (e6)~ 2 ] or superdeformed [~5 MeV _1 (e6)~ 2 ] rotational bands. 

Similar as ferromagnetism and antiferromagnetism in condensed matter physics, the ex- 
istence of the antimagnetic rotation (AMR) in nuclear physics is also an interesting issue. 
For an antiferromagnet, one-half of the atomic dipole moments are aligned on one sublattice 
and the other half are aligned in the opposite direction on the other sublattice. In such 
way, the net magnetic moment in an antiferromagnet is canceled out. However, it is still an 
ordered state since the isotropy of such a state is also broken 

In analogy with an antiferromagnet, antimagnetic rotation 
some specific nearly spherical nuclei, in which the subsystems of valence protons (neutrons) 
are aligned back to back in opposite directions and nearly perpendicular to the orientation 
of the total spin of the valence neutrons (protons). Such arrangement of the proton and 
neutron angular momenta also breaks the rotational symmetry in these nearly spherical 
nuclei and causes excitations with rotational character on top of this bandhead. Along this 
band, energy and angular momentum are increased by simultaneous closing of the two proton 
(neutron) blades toward the neutron (proton) angular momentum vector. Consequently, a 
new kind of rotational bands in nearly spherical nuclei is formed showing some analogy with 
an antiferromagnet. A schematic illustration for AMR is shown in Fig. 

The experimental indicators for antimagnetic rotation can be summarized as: 

1) a AI = 2 sequence with E2 transitions only, as the cancellation of the magnetic moments 
leads to the absence of the Ml transition; 

2) weak E2 transitions corresponding to the small deformation parameter < 0.15, and 
the B(E2) values decreasing with spin; 

3) a smooth variation in the 7 transition energy with angular momentum; 

4) a substantial dynamic moment of inertia, similar to the MR band, corresponding to the 
large ratio of the jW/B(E2) > lOOMeV" 1 (eb)~ 2 . 

The MR and AMR 8Hl0j discussed above are based on the assumption that the nucleus 
involved is not triaxially deformed. If the nucleus is triaxially deformed and the aplanar 



rotation is then allowed, another exotic phenomenon, nuclear chirality, may occur 



15 



16|. 
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FIG. 2: A schematic illustration for the spin-coupling scheme of antimagnetic rotation. Instead of 
the shears mechanism for magnetic rotation, the two shearslike configurations here result that the 
magnetic moments are antialigned and cancel, and therefore the B(M\) values vanish. 



B. Experimental progress 



Following the observation of many long cascades of mag netic dipole 7-ray transitions in 



17H23l|. several attempts 



24- 



have 



the neutron deficient Pb nuclei in the early 1990s 
been made to measure the lifetimes for the states in MR bands. The high-accuracy lifetime 
measurements for four Ml bands in 198 > 199 Pb performed with GAMMASPHERE provided 



a clear evidence for shears mechanism (29]. Additional evidence for t 



le shears mechanism 



was provided by measuring the ^-factor of a dipole band in 193 Pb 30|. It is demonstrated 



that at the bandhead the longitudinal component of the magnetic moment n\\ has the value 
expected for an opening angle of 90° of the two blades composed of the suggested particles 
and holes. 

From then on, more and more MR bands have been observed not only in the mass region 
of A ~ 190 but also in A ~ 80, A ~ 110, and A ~ 140 regions. A compilation in these 
four mass regions up to December 2006 including 178 magnetic di pol e bands in 76 nuclides 



is given in Ref. 31] and additional data can be found in Refs. 32 



4311 . Recent observations 



in 58 Fe 



44J and 60 Ni 45j have extended the observed MR mass region to A = 60. In total, 



more than 195 magnetic dipole bands spread over 85 nuclides have been observed, which 
have been summarized in the nuclear chart in Fig. [3j 
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It is difficult to define how much contribution from the collective rotation can be expected 
in the MR bands. Obviously either the alignment of the angular momenta for the valent 
nucleons or the collective rotation due to deformation will cost energy. For a given angular 
momentum, the nuclear system will try to minimize its energy via a competition between 
the kinetic energy due to the collective rotation and the potential energy due to the closing 
of the blades of the shears blades. Weaker quadrupole deformation will lead to purer MR 
band and larger quadrupole deformation will reduce the B(M1)/B(E2) ratios. Therefore 
it might be possible to observe the competition and transition between the electric and the 
magnetic rotations. The experimental information for the lifetime measurement in 106 Sn has 



suggested an extremely large ratio / B{E2) > lOOOMeV 1 (eb] 
example of almost pure magnetic rotation in a spherical nucleus 



"\ which may provide an 
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FIG. 3: The nuclides with magnetic rotation observed in the nuclear chart. The red squares 
represented 56 nuclides with 120 magnetic dipole bands observed before 2000 and compiled in 
Ref. 



31] . The green squares represented the corresponding data after 2000. 

AMR is expected to be observed in the same regions as MR in the nuclear chart (9). How- 
ever, it differs from magnetic rotation in two aspects. Firstly, there are no Ml transitions in 
the AMR band since the transverse magnetic moments of the two subsystems are antialigned 
and canceled out. Secondly, as the antimagnetic rotor is symmetric with respect to a rota- 



6 



tion by it about the rotating axis, the energy levels in AMR band differ in spin by 2h and are 
connected by weak E2 transitions reflecting the nearly spherical core. Moreover, the AMR 
phenomenon is characterized by a decrease of the B(E2) values with spin, which has been 



confirmed by lifetime measurements 



471 ] . To date, AMR has attracted lots of attentions and 



has been observed in Cd isotopes including Cd [48|], Cd 



no 



Cd 



108 



Cd 



49 



50| i° 9 Cd |5l|, 



52|], and 112 Cd [53j. The other candidates include 100 Pd [54j, and 144 Dy [55]. 



C. Theoretical progress 

From the theoretical side, both the MR and AMR bands have been extensively discussed 
in recent years. Using a semiclassical particle plus rotor model based on angular momentum 



geometry 
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56 



571 ]. the competition between shears mechanism and core rotation has 
been investigated. As the mean field approach is easy to construct classical vector diagrams 
showing the angular momentum composition, it is widely used to understand the structure 
of these rotational bands. Since for the magnetic rotation bands, the axis of the uniform 
rotation does not coincide with any principal axis of the density distribution, a description 
of these bands requires a model going beyond the principal axis cranking, which leads to 
the development of the tilted axis cranking (TAC) approach. 

The semi-classical mean field approximation for tilted axis rotation can be traced back 
to the 1980s 



OS 



59j. After the first self-consistent TAC solutions were found in Ref. [11] . 



the qualities and interpretations of the TAC approximation were discussed and examined in 
Ref. |60J with the particle rotor model. Based on the TAC approximation, lots of applications 
are carried out in the framework of the pairing plus quadrupole model or the shell correction 
method ^Q- 

In the above mentioned TAC investigations, however, the polarization effects which are 
expected to strongly influence the quadrupole moments and thus the B(E2) values, are 
either completely neglected or taken into account only partially by minimizing the rotating 
energy surface with respect to a few deformation parameters. Moreover, the nuclear currents, 
which are the origin of symmetry violation in nuclei with MR and/or AMR, are not treated 
in a self-consistent way in these models. Therefore, it is necessary to investigate the MR 
and AMR based on thoeries which describe the polarization effects and nuclear currents 
self-consistently without additional parameters. Such calculations are not simple, but they 
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are nowadays feasible in the framework of density functional theory (DFT). 

The DFT with a small number of parameters allows a very successful description of 
ground-state and excited-state properties of nuclei all over the nuclear chart. On the basis 
of the density functionals, the rotational excitations have been described in the principal- 



axis cranking Hartree-Fock 
zero-range Skyrme force 



62 



TF) or Hartree-Fock-Bogoliubov (HFB) framework with the 



64J as well as the density dependent Gogny force [65] non- 



relativistically. The three-dimensional TAC version, which allows to study nuclear aplanar 
rotation, has been developed with the Skyrme density functional and applied to not only the 



magnetic rotation in Gd in Ref. 



66] but also chirality (15J in A ~ 130 mass regions 
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The covariant version of DFT takes the Lorentz symmetry into account in a self-consistent 
way and has received wide attention due to its successful description for lots of nuclear phe- 
nomena in stable as well as exotic nuclei 69M7l| . The representations with large scalar and 
vector fields in nuclei, of the order of a few hundred MeV, provide more efficient descriptions 
than non-relativistic approaches that hide these scales. It can obtain reasonable nuclear sat- 
uration properties in infinite nuclear matter with Brueckner method 72|, [73] , reproduce well 
the measurements of the isotopic shifts in the Pb region 74] . give naturally the spin-orbit 



potential, explain the origin of the pseudospin symmetry 
try 



75 



77 



83] and the spin symmetry in the anti-nucleon spectrum 



76 1 as a relativistic symme- 



84 



and is reliable for 



nuclei far away from the ^-stability line 



N6 



87], etc. Moreover, it is of particular im por tance 



that the covariant density functional theory (CDFT) includes nuclear magnetism 88], and 
provides a consistent description of currents and time-odd fields, which plays an important 
role in the nuclear rotations. 

Based on the CDFT, Koepf and Ring 



891 ] first developed the principal axis cranking 



relativistic Hartree calculations. Konig and Ring [90] included the effects of isoscalar and 
isovector baryon currents that are the sources of the magnetic potentials in the Dirac equa- 
tion. The nuclear magnetism appeared to be crucial for the quantitative understanding of 
nuclear moment of inertia. Large scale cranking relativistic Hartree-Bogoliubov calculations 
using the central part of the Gogny force with the PIS parameter set acting in the pairing 



channel has been performed and summarized in Ref. [91] . These calculations went beyond 
the mean-field by invoking an approximate number projection which appears to have great 
impact on the values. 

The three-dimensional cranking CDFT has been developed in Ref. 92], in which the 



S 



nucleon and meson fields are expanded in terms of three-dimensional harmonic-oscillator 
eigenfunctions. In comparison with the principal axis cranking (PAC) case, the signature 
symmetry is broken and the parity is thus the only symmetry left 92j. Even without the 
pairing correlations, it is very time consuming to perform a three-dimensional cranking 
calculation. Therefore the three-dimensional cranking CDFT is firstly used to examine its 
applicability to the shears bands, i.e., the magnetic rotation in 84 Rb 921 ] . 

Focusing on the magnetic rotational bands, a completely new computer code for the 
self-consistent two-dimensional cranking CDFT based on the non-linear meson-exchange 



interaction has been established 



931 ] . This new code includes significant improvements, such 



as implanting the simplex symmetry V y T, quantum number transformation, and orientation 
constraints. Compared with the three-dimensional cranking CDFT, the computing time for 
this two-dimensional cranking CDFT model has been considerably reduced and thus it allows 
systematic investigations and has been applied to heavy nucleus like 142 Gd 93]. 

The CDFT with the point-coupling interaction become very popular in recent years 
owing to the following advantages, 1) it avoids the possible physical constrains introduced 
by explicit usage of the Klein-Gordon equation, especially the fictitious a meson; 2) it is 



possible to study the naturalness of the interaction 



94 



95] in effective theories for related 



nuclear structure problems; 3) it is relatively easy to include the Fock terms 96| and to 



investigate its relationship to the nonrelativistic approaches 



97|. 



The TAC model based on the CDFT with the point-coupling interaction was established 



in Ref. 



98 



which introduces further simplification and reduces computing time. It has been 



applied successfully for the MR ran ging from light nuclei such as Ni 



and 58 Fe 



medium heavy nuclei such as In |40j], and to heavy nuclei such as 



198,199 



Pb 



to 



9J|. It also 



provides a fully self-consistent and microscopic investigation for the AMR in 105 Cd 



More examples of its application for the AMR bands in 112 In can be found in Ref. 53] . 



100]. 



In the review, we will mainly focus on the microscopic CDFT investigation for the MR 
and AMR. In Section II, the theoretical framework of the tilted axis cranking CDFT with 
the point-coupling interaction will be introduced. In Section III, taking the seminal MR 
in 198 Pb as an example, the CDFT description for the MR including the single particle 
Routhian, energy spectra, deformation evolution, shears mechanism, electric and magnetic 
transitions properties, etc. will be outlined. Similar applications for MR in the mass regions 
A ~ 60, 80, 100, 140 are summarized in the next Section. In Section V, taking 105 Cd as an 
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example, the CDFT investigation for the AMR are discussed. Finally, the summary and 
perspectives are given in Section VI. 



II. TILTED AXIS CRANKING CO VARIANT DENSITY FUNCTIONAL THE- 
ORY 



The basic idea of the cranking model is based on the following classical assumption |101| : 
If one introduces a coordinate system which rotates with constant angular velocity around 
a fixed axis in space, the motion of the nucleons in the rotating frame is rather simple if 
the angular frequency is properly chosen. In particular, the nucleons can be thought of as 
independent particles moving in an average potential which is rotating with the coordinate 
frame. 

The advantages of the cranking model include, 1) providing a microscopic description for 
rotating nucleus; 2) describing the collective angular momentum as the sum of the single 
particle ones; 3) working extremely well at high angular momentum where the assumption 
of uniform rotation applies. 

For the magnetic rotations, the axis of the uniform rotation does not coincide with any 
principal axis of the density distribution, a description of the AI = 1 rotational bands 
requires going beyond the principal axis cranking, which leads to the TAC approach 11]. 
The initial version of the TAC approach was developed based on the potentials of the Nilsson 
type, which are combined with a pairing plus quadrupole model or the shell correction 
method for finding the deformation. 

In this section, we first present the basic idea of TAC approach by taking the simple 
pairing plus quadrupole model fPQTAC) as an example. The details of the PQTAC model 



can be found in Refs. 



6lfl . In Subsection III Bl the theoretical framework of the 
covariant density functional theory will be briefly introduced. The combination of the TAC 
approach and the covariant density functional theory, i.e., the so-called tilted axis cranking 
relativistic mean-field theory, is given in the Subsection III CI Finally, Subsection III Dl will be 
devoted to the numerical techniques. 
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A. Tilted axis cranking model 



Assuming that a nucleus rotates around an arbitrary axis with a uniform velocity f2, the 
corresponding Hamiltonian (with the pairing plus quadrupole interaction) in intrinsic frame, 
i.e., the Routhian, can be written as, 

2 

H' = H-n-j = H sph -^Yl Q\Q» - Gp]p - a# - n • j. (i) 

p,=-2 

The constraint Q ■ J ensures that the rotational states have a finite angular momentum J. 
The justific atio n of pairing plus quadrupole Hamiltonian is described in many textbooks 



(e.g., Ref. [101| ). In PQTAC model, H sp h is approximated by the spherical part of the 
standard Nilsson Hamiltonian. The residual interaction includes two parts. The quadrupole 
interaction defined by the quadrupole operators is responsible for the quadrupole deforma- 
tion of the mean field. The pairing interaction takes a monopole form. The term XN controls 
the particle number N. The Hamiltonian ([1]) is written only for one kind of nucleons. The 
terms H sp h and should be understood as sums of a proton and a neutron part, and there 
are terms —GP^P and — XN for both protons and neutrons. In practice, the actual values 
of the force constants x an d G depend on the configuration space under consideration and 
are determined from experimental data. 

Approximating the nuclear wave function as a Slater determinant in the quasiparticle 
space and neglecting the exchange terms, the single particle Routhian h! becomes 

i 2 

h' = h sph --J2 (g M Qi + g*g At )-A(pt + p)-AiV-n. J (2) 

fi=-2 

with the deformed potential g M = x(Qn) an d the pairing potential A = G(P). The Fermi 
surface A is determined by N = (N), with N the particle number of neutrons or protons. 

For given parameters and A, solving Eq. (J2J) is similar to solving a cranking Nilsson 
Hamiltonian. Here and A depend on the wave function and are determined by itera- 
tion, or equivalently by minimizing the total Routhian E'^q^, A) = (H') for each rotational 



frequency Q. It has been shown in Ref. [58| that a self-consistent solution is equivalent to 
the requirement that the vector of the rotational frequency f2 is parallel to the vector of 
the angular momentum. Finally, the total energy as function of the angular momentum is 
obtained by 

E(I) = E'(Ct) +QI(n), (3) 
11 



where is determined with I = \{J)\. 

Usually, it is convenient to choose the principal axes of the density distribution as the 
axes of the coordinate frame, i.e., 

9-1 = ?i = 0, q- 2 = <?2- (4) 
In this case, the deformation of the potential is specified by two intrinsic quadrupole moments 



gn and qi. Therefore, it is useful to reformulate the TAC model in this intrinsic frame 
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9. 



Cranking model are based on the classical treatment of the total angular momentum and 
the assumption of uniform rotation, which have the consequence that angular momentum 
conservation is violated. The connection with the quantal spectra is made by means of semi- 
classical expressions for the energy and transition matrix elements. Hence, it is important 
to investigate how well these approximations work for the description of the experimental 
observables. Moreover, in the TAC approach, the rotational axis does not coincides with any 
principal axis, i.e., it is tilted away from the principal axis. Different from the principal axis 
cranking model, the signature is no longer a good quantum number in the TAC. This leads 
to a problem of how to interpret the TAC solutions and construct the excitation spectrum 
from the TAC quasiparticle levels avoiding spurious states. 

The above questions have been discussed in Ref. 6(| by comparing with the particle rotor 
model, which treats the quantal angular momentum dynamics properly. It was found that 
the TAC approach quantitatively accounts both for the energies and the intra band transition 
rates of the lowest bands generated by one or two quasi particles coupled to an axial rotor. 
The TAC provides an accurate description of the bandhead, except in cases, when substantial 
alignment of quasi particle angular momentum occurs at very low frequency. 

The TAC approach is based on the mean-field theory. Therefore, one may easily study 
multi-quasiparticle excitations, and the consequences of changes of the deformation or the 
pairing. It gives transparent classical vector diagrams pictures of the angular momentum 
coupling, which is of great help to understand the structure of rotation bands. The drawback 
of the TAC model is that it cannot describe the gradual onset of signature splitting as well 
as the mixing of bands with substantially different quasiparticle angular momentum as in 
the standard cranking theory. 
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B. Covariant density functional theory 



The covariant density functional theory can be traced bac k to t he s uccessful relativistic 
mean- field (RMF) models introduced by Walecka and Serot 102 . Il03 ], which was further 
developed and applied by many groups [69H7l|. The most popular RMF models are based 
on the finite-range meson-exchange representation, in which the nucleus is described as a 
system of Dirac nucleons that interact with each other via the exchange of mesons. The 
nucleons and the mesons are described by the Dirac equation and the Klein-Gordon equation, 
respectively. The isoscalar-scalar a meson, the isoscalar-vector u meson, and the isovector- 
vector p meson build the minimal set of meson fields that, together with the electromagnetic 
field, is necessary for a description of bulk and single-particle nuclear properties. Moreover, 
a quantitative treatment of nuclear matter and finite nuclei needs a medium dependence 
of effective mean-field interactions, which can be introduced by including nonlinear meson 
self-interaction terms in the Lagrangian or by assuming explicit density dependence for the 
meson-nucleon couplings. The detailed formulism of the meson-exchange representation of 



CDFT can be found in Refs. 



3- Ill- 



More recently, this framework has been reinterpreted by the relativistic KohnCSham 
density functional theory, and t he functio nals have been developed based on the zero-range 



point-coupling interaction 



*6 



104 



105], in which the meson exchange in each channel 



(scalar-isoscalar, vector-isoscalar, scalar-isovector, and vector-isovector) is replaced by the 
corresponding local four-point (contact) interaction between nucleons. In recent years, the 
point-coupling model has attracted more and more attention owing to the following ad- 
vantages. First, it avoids the possible physical constrains introduced by explicit usage of 
the Klein-Gordon equation to describe mean meson fields, especially the fictitious a me- 
son. Second, it is possible to study the role of naturalness 94j, |95( in effective theories for 



96|, 



nuclear-structure-related problems. Third, it is relatively easy to include the Fock terms 
and provides more opportunities to investigate its relationship to the nonrelativistic ap- 
proaches 97J . In the following, we present the theoretical framework of the CDFT with the 
point-coupling interaction. 

The basic building blocks of the covariant density functional theory with point-coupling 
vertices are 

(i/jOTi/j), OE{l,r}, r6{l, 7p) 7 5 ,7 5 7„ )V }, (5) 
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where ip is Dirac spinor field of nucleoli, f is the isospin Pauli matrix, and V generally 
denotes the 4x4 Dirac matrices. There are ten such building blocks characterized by their 
transformation characteristics in isospin and Minkowski space. In the following, the vectors 
in the isospin space are denoted by arrows and the space vectors by bold type. Greek indices 
fi and v run over the Minkowski indices 0, 1, 2, and 3. 

A general effective Lagrangian can be written as a power series in ipOTip and their 
derivatives, with higher-order terms representing in-medium many-body correlations. In 
the actual application we start with the following Lagrangian density of the form: 

£ = C iTCC + £ 4f + £ hot + C dcr + £ cm 

= ip(i'-f fl d^ — m)tp 

1 - - 1 1 - - 1 - 

—Syd^^^Y^) - \5 TS d v {^)d v {^) - ^Tvd^^^^r^) 
--F^F» V - el^pL^A^ (6) 

which includes the Lagrangian density for free nucleons £ free , the four-fermion point-coupling 
terms £ 4f , the higher order terms £ hot accounting for the medium effects, the derivative terms 
£ der to simulate the effects of finite-range which are crucial for a quantitative description for 
nuclear density distributions (e.g., nuclear radii), and the electromagnetic interaction terms 
£™. The higher order terms lead in the mean field approximation to density dependent 
coupling constants with a density dependence of polynomial form. 

For the Lagrangian density in Eq. , m is the nucleon mass and e is the charge unit for 
protons. and are respectively the four-vector potential and field strength tensor of 
the electromagnetic field. There are totally 11 coupling constants, as, ay, a rs, oltv-, As, 7s, 
7v, 5s, <5y, Sts, an d Stv, m which a refers to the four-fermion term, (3 and 7 respectively the 
third- and fourth-order terms, and S the derivative couplings. The subscripts S, V, and T 
respectively indicate the symmetries of the couplings, i.e., S stands for scalar, V for vector, 
and T for isovector. The isovector-scalar channel including the terms cxts and Sts m Eq. (JH]) 
is usually neglected since a fit including the isove ctor- scalar interaction does not improve 



the description of nuclear ground-state properties 1051 ] . Furthermore, the pseudoscalar 75 



and pseudovector 757^ channels are also neglected in Eq. §6§ since they do not contribute 
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at the Hartree level due to parity conservation in nuclei. 

The Hamiltonian density, i.e., the 00 components of the energy- momentum tensor can be 
obtained by the Legendre transformation 

% = T 00 = ^ - £, (7) 
where 0, represents the nucleon or photon field. Therefore, the total Hamiltonian reads 



H 



d 3 xU 



J d 3 x {i> [-ij ■ V 



m] ip 



+^as(#)(#) + ^«y(^7^)(^7^) + ]^ TV {ipf^^){ipf^i)) 



-5 TV [d (iJTl^)d°(H^) + V(^tV/0 • V^tY 



One can describe the nucleon field as 



(8) 



k 



(9a) 
(9b) 



where ipk{x) and ip' k (x) construct a complete set of the Dirac spinor. The operators at and 
a k respectively represent the creation and annihilate operators for nucleons, whereas d k and 
d\ represent those for the antinucleons which corresponds to the negative energy states in 
the Dirac sea. Usually the contributions from the Di rac s ea are neglected in the framework 
of CDFT, i.e., the so-called "no-sea" approximation 102 ]. 



The nuclear ground-state wavefunction |$) is assumed as a Slater determinant 



i*>=iI fl H->' 



(10) 



k=l 
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where |— ) is the physical vacuum. This leads to the replacement of the operators ^(OT)^ 
in Eq. fl6]) by their expectation values which become bilinear forms of the nucleon Dirac 
spinor ip/, in the Hartree approximation, 

i>(6T)^ -> ($| : ^{OV)^ : |$) = Y M^(OT)iM^), (n) 

k>0 

where i indicates S, V, and TV, and the sum runs over only positive energy states. 
Accordingly, the energy density functional for a nuclear system can be represented as, 

Eqdf = (&\H\$) 

/[ A 11 1 -> 

d 3 x < Y ii ( a ' P + P 171 ) + 2«5Pl + -jxvitif + i^tv6tv)Jtv 
I k=l 

+^SP 3 S + \lsPs + \lv{3»f? - \h [d psd°p s + Vp 5 ■ Vp 5 ] 

-\&v [doJ,d°f + Vj, ■ Vf] - l -5 TV [d (jTv)»d ft v + V(j TV ), ■ Vj£ y 

+eApj$ - F^d°A, + \f, u F^ , (12) 

where the local densities and currents read, 



Ps{x) 



($| : W : |$) = Mx)Mx), (13a) 



k=i 

A 



f{x) = ($| : fy*V : 1$) = ^^(^7^^), (13b) 

fc=i 

A 

j£ v (x) = ($| : ■ 1$) = ^^(x) 7 M r^(x), (13c) 

k=l 
k=l 

For the stationary case, ipk{x) = ipk{ r )e~ lSkt , by minimizing the energy density functional 
Eq. ffl2|) with respect to tpk, one obtains the Dirac equation for the single nucleons 

[a • (-iV - V) + V + /3(m + = (14) 

Here, the local scalar S{r) and vector V^{r) potentials read 

S(r) = E s , V^r) = W + r-^ TV , (15) 



1(3 



where the nucleon scalar-isoscalar Eg, vector-isoscalar E M , and vector-isovector E^ self- 
energies are given in terms of the various densities, 

£s = asps + PsPs + ISPs + 5 s &ps, (16a) 
^ = av f + 7 y(j M ) 3 + S v Af + e^^A", (16b) 

Z>TV = a TV3TV + StvAJtv ( 16C ) 

Similarly, one can also obtain the Coulomb field which is determined by Poisson's equa- 
tion 

- M"(r) = ef p {r). (17) 

For a system with time reversal invariance, the space-like components of the currents ji 
in Eq. ffl3|) and the vector potential V(r) in Eq. ffl5|) vanish. Further assuming that the 
nucleon single-particle states do not mix isospin, i.e., the single-particle states are eigenstates 
of r 3 , only the third component of isovector potentials E^, y survives. The total energy is 
thus given by 



E, 



CDF 



/[ - 111 
d 3 x < ^2 ( a ■ P + P m ) + ^ a sP 2 s + i^vPv + ^tvPtv 
I k=i 



+^PsP s + ^7sP| + \lVPv + ^3p S Ap S 

+^5 v p v Ap v +^5 T vpTvAp T v + ^eAoPpj . (18) 

As the translational symmetry is broken in the mean-field approximation, the center-of- 
mass (cm.) correction should be made for the spurious cm. m otion. N owadays, this is 



usually done by including the microscopic cm. correction energy 106l . Il07j ] 



with A being the mass number and P c . m . = Yli pi being the total momentum in the cm. 
frame. Therefore, the total energy for the nuclear system becomes 

E to t = E C df + E™£. (20) 

In the framework of the CDFT with the point-coupling interaction, a new point-coupling 
effective interaction PC-PK1 has been proposed by fitting to observables of 60 selected spher- 



ical nuclei, including the binding energies, charge radii, and empirical pairing gaps 
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86|]. This 



effective interaction particularly improves the description for isospin dependence of binding 
energies and it has bee n su ccessfully used in des cribing the Coulomb displacement energies 
between mirror nuclei 



108 



fission barriers 



109], the new and accurate mass measurement 



results at Gesellschaft fiir Schwerionenforschung (GSI) 87|, etc. 

By using the PC-PK1, the masses of nuclei with neutron n umbe r N > 8 and proton 
number Z > 8 in the atomic mass evaluation of 2003 (AME03) llCll ] have been calculated 
in the framework CDFT with the axial symmetry. The mass differences between the exper- 
imental data and the calculated results are shown in Fig. HI It is found that the CDFT with 
PC-PK1 can reproduce the experimental data quite well and the corresponding root-mean- 
square (rms) deviation is 1.422 MeV, which is much smaller than therms deviation value of 



2.25 MeV given by the meson-exchange effective interaction TMA [111] - Note that in the re- 
sults given by PC-PK1, the rotational correction energy due to the violation of the rotational 
symmetry is considered for deformed nuclei with the moment of inertia X > 3.5 h 2 /MeV 
by E rot = —^z(J 2 ), where J is the angular mome ntum ope rator and X is the moment of 
inertia calculated from the Inglis-Belyaev formula 112l - lll4j . The parameter b is chosen 
as 1, 1.4, and 1.6 for even-even, odd-A, and odd-odd nuclei, respectively. As discussed in 



Ref. 



87| . further improvements on a proper treatment of the rotational correction energy 



for odd nuclei is necessary. 
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FIG. 4: Mass differences between the experimental data [l lQl ] and the covariant density functional 
theory calculations with the point-coupling interaction PC-PK1 [8Q]. 
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C. Tilted axis cranking CDFT 



In order to describe the magnetic rotations in a microscopical and self-consistent way 
the tilted axis cranking CDFT should be developed. This approach has been realized based 
on either the meson-exchange interaction 



92 



93] or the point-coupling interaction [98 ]. 



The TAC-RMF approach based on the meson-exchange interaction is described in detail in 



Ref. 



93] . In the following, the main discussion will be focused on the TAC-RMF approach 



based on the point-coupling interaction 



98 



Assuming that the nucleus rotates around an axis in the xz plane, the Lagrangian in 
Eq. is transformed into a frame rotating uniformly with a constant rotational frequency, 



ft = (Q x ,0,Q z ) = (ncos0n,O,fisin0 r 



(21) 



where 8q := <(0, e x ) is the tilted angle between the cranking axis and the x axis. From this 
rotating Lagrangian, the equation of motion for the nucleons can be derived equivalently by 
either starting from a special relativistic transformation 88], or more generally adopting the 
tetrad formalism in the framework of general relativistic theory 1 151 ] . One can thus finds 

[a • (-iV - V) + p{m + S) + V - n ■ Jtyk = e^A, (22) 

where J = L + is the total angular momentum of the nucleon spinors, and represents 
the single-particle Routhians for nucleons. The relativistic fields S(r) and V^{r) read 

S(r) = a s p s + PsP 2 s + lsP 3 s + $s&Ps, (23a) 
V(r) = a v p v + 'yvpy + 5 v Ap v + T 3 a TV p TV + t 3 5tv^Ptv + eA°, (23b) 
V(r) = a v j v + -f V (j v ) 3 + 5 v Aj v + r 3 a T v3Tv + T 3 5 T vAjTv + eA. (23c) 

As usual, it is assumed that the nucleon single-particle states do not mix isospin, i.e., the 
single-particle states are eigenstates of t$. Therefore only the third component of isovector 
potentials survives. The Coulomb field A°(r) is determined by Poisson's equation 



AA°(r 



ep c 



(24) 



The spatial components of the electromagnetic vector potential A(r) are neglected since 



their contributions are extremely small 88 



Since the Coriolis term O • J in the Dirac equation (122]) breaks time reversal symmetry 
in the intrinsic frame, the nucleon currents are induced and as a consequence the spatial 
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components of the vector potential V(r). The densities and currents in Eqs. ( |23|) have the 
form 

A 



Ps\T 



i=l 
A 

Pv(r) = 



i=l 
A 



Jv{r) 



i=l 
i=l 



pTv{r) 



i=l 
A 



P 



(r) = $>1 



r)aV'i(T-), 

r)r 3 ^i(r), 

r)m- 3 V>i(r), 
a-r 3 



(25a) 
(25b) 
(25c) 
(25d) 
(25e) 
(25f) 



i=l 



Here, the "no-sea" approximation is also adopted, i.e., the sums run over only the particles 
states in the Fermi sea and the contribution of the negative-energy states are neglected. 

By solving the equation of motion iteratively, one finally obtains the total energy in the 
laboratory frame 

-Etot = E kin + E int + E cou + E c , m _, (26) 
which is composed of a kinetic part 

r A 

E kin = d 3 rJ2 " P + (27) 



an interaction part 

E- mt 



i=l 



d 3 r S^a s p 2 s + ^f3 s p 3 s + \lsPs + ^sPs^Ps 

+ ^ a v(pv ~ 3 ■ j) + 7}Utv(Ptv ~ 3tv ■ hv) 
+ \lv(p 2 v ~ 3 ■ 3? + -&v{pv&Pv - jAj) 



+ -5tv(ptv<Aptv - jrvAjrv] 



an electromagnetic part 



E, 



d r-eA p c , 



(28) 



(29) 
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and the center-of-mass (cm.) correction energy E c . m . accounting for the treatment of center- 
of-mass motion. 

For each rotational frequency Q, the expectation values of the angular momentum com- 
ponents J = (J x ,J y ,Jz) i n the intrinsic frame are given by 

A 

J, = (4> = (30a) 

8=1 

J y = 0, (30b) 

A 

J z = (l) = Y,^\ (30c) 
i=i 

and by means of the semiclassical cranking condition 

J= y/ A) 2 + (I) 2 = v^+l), (31) 

one can relate the rotational frequency Q to the angular momentum quantum number / in 
the rotational band. 

The orientation of the angular momentum vector J is represented by the angle 8j : = 
<( J, e x ) between the angular momentum vector J and the x axis. As mentioned before, in 
a fully self-consistent calculation, the orientation 9j of the angular momentum J should be 
identical to the orientation 6q of the angular velocity ft. 

The quadrupole moments Q20 and Q22 are calculated by 



k/i!^ 3 * 2 - r2 >> ( 32a ) 



Q22 = v^ 2 -^' ( 32b ) 

and the deformation parameters (3 and 7 can thus be extracted from 

P = \l a io + 2 «22 5 (33a) 
7 = arctan 



O20 



(33b) 



by using the relations 



3A 

Q20 = ^r^o a 20, (34a) 
3A 

Q22 = ^^0«22, (34b) 
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with R = 1.2 A 1 / 3 fm. Note that the sign convention in Ref. [10 1| is adopted for the 
definition of 7 here. 

The nuclear magnetic moment, in units of the nuclear magneton, is given by 



i=i 



mc 
he 



qi/jj(r)r x aipi(r) + fc0j(r)/3S^)j(r) 



(35) 



where the charge q (q p = 1 for protons and q n = for neutrons) is in units of e, m the 
nucleon mass, and k the free anomalous gyromagnetic ratio of the nucleon (n p = 1.793 and 
K n = -1.913). 

From the quadrupole moments and the magnetic moment, the B(M1) and B(E2) tran- 
sition probabilities can be derived in semiclassical approximation 



B(M1) = — fi ± = — (fi x sin 9 j — ji z cos^j 



57T 



57T 



B(E2) 



8 



g^ cos 2 ^ + A /-g^ 2 (l + sin 2 ^) 



(36a) 
(36b) 



where an( ^ Q22 corresponds to the quadrupole moments of protons. 



D. Numerical techniques 

Orientation constraint In the usual PAC programs (one-dimensional cranking), the 
principal axes of the densities and fields are implemented to be along the x, y, and z axis. 
For the TAC code (two-dimensional cranking), it allows for arbitrary rotations of the density 
distributions around the intrinsic y axis. The freedom of rotations around the y axis can 
lead (in particular for Q = and for small Q values) to instabilities during the iterative 
solution because the solutions with different orientations in the xz plane are degenerate. 
Therefore, the x, y, and z axes are enforced to be identic al w ith the principal axes of the 



density distribution by introducing a quadratic constraint 101] for the expectation value of 
the quadrupole moment 



i.e., by minimizing 



(Q 2 - 1 ) = -\^(xz)=0, (37) 

V o7T 



{H') = {H)+ l -C{{Q^ 1 )-a^ 1 )\ (38) 



with d2-i = 0, and C being a spring constant, which, if properly chosen, has no influence 
on the final result. 
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Expansion in harmonic oscillator basis In the code for the solution of the relativistic 
tilted axis cranking equations, the Dirac spinors are expanded in terms of three-dimensional 
harmonic oscillator wave functions in Cartesian coordinates, 



\n x ,n y ,n z ) = (p nx (x)<p ny (y)(p nz {z). 



(39) 



The normalized oscillator function tp nk (xk) in fc-direction {x^ = x,y,z) are given by 



N„, ,%k\ r ^-f X k\2l 

fe v— )exp[--(— ) ], 



VOfc °k * Ok 



(40) 



where N n = (y/n2 n n\) 1//2 is a normalization factor and 



2 d n _c2 



■1 )"e< — e 



(41) 



are the Hermite polynomials 116]. 

We can take the following basis states 



<p a [T,s) 



(r,s\a) = i nv ¥n x { x )¥n y (y)VnM 



V2 



1 

_l\n x +l 



(42) 



and 



<Pa(r,s) = (r,s\a) = {-i) ny <Pn x { x )Vn y {y)<PnA z ) 




(43) 



which correspond to the eigenfunctions of the simplex operation with the positive and 
negative (— i) eigenvalues, respectively. The pha se fa ctor i ny has been added in order to 



have real matrix elements for the Dirac equation 
The Dirac spinor for the nucleon has the form 



Xi(t), 



(44) 



where Xi{t) is the isospin part. In the tilted axis cranking calculations, the simplex symmetry 
is violated, and therefore the large and small components of the wave function in Eq. ( 144)) 
have to be written as linear combinations of the sets fT42l and (T4"3"l) with different simplex: 



a a 



(45) 
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Since the large and small components in the Dirac equation have different parity, the sums in 
the expansions for the large and the small components have to run over oscillator quantum 
numbers with even N = n x + n y + n z or odd N respectively. This is indicated in Eq. (113]) 
by the indices a and a. 

On this basis, the solution of Dirac equation ([22]) is obtained by the matrix diagonalization 



(46) 





( fa\ 




( fai\ 


n 


fat 




fat 




= Si 




9ai 




9ai 




\9m J 




\93i J 



where the Hamiltonian matrix H has the form 



/ 



(a\M* + V - nj\a') {a\M* + V -nJ\a') (o|<r(V - iV)\a!) 



(a\a(V -iV) 
(a\a{V -iV) 



a 



a 



(a\M* + V - QJ\a') (a\M* + V - £U\a') (a|<r(V - iV) \a!) 

(d|-cr(V -iV)\a') (5|-<r(V - iV)\a') (a\-M* + V - QJ\a') (a\-M* + V - ftJ 
\ <q|-ct(V -iV)\a!) (a\-o-(V - iV)\a') (a\-M* + V - OJ|a') (a\-M* + V -Clj 



I 

(47) 



Note that the Coriolis term OJ breaks the invariance with respect to time reversal and 
with respect to rotations around the well as around the z axis. Therefore, only 

the invariance of space reflection V and the combination of time reversal and reflection in y 
direction V y T are valid and used in the code. 

For the evaluation of the Coulomb field, due to its long range character, an expansion 
in harm onic oscillator states is very difficult and therefore the standard Green function 



method 



1 18l | is used for the calculation of the Coulomb field in each step of the iteration. 



Configuration constraint Normally the rotation bands are built on specific proton and 
neutron configurations. The orbits to be blocked are usually given in the spherical basis, 
i.e., by the spherical quantum numbers \nljm). The equations of motion are solved by 
expanding the Dirac spinors in terms of the three-dimensional harmonic oscillator functions 
in the Cartesian basis Eq. (T45]) labeled by the quantum number \n x ,n y ,n z ,±i). In order 
to identify which orbits have to be blocked, one need to transform the wave functions from 
the Cartesian basis with the quantum number \n x n y n z ± i) to a sp h erica l basis with the 
quantum numbers \nljm) using the techniques given in Refs. 119l . Il20| . Consequently, 
one can block the levels which have the maximal overlap with the required \nljm) orbits. 
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These techniques are considerably simplified, if we work in an isotropic Cartesian basis with 
identical basis parameters b x = b y = b z = b in Eq. (j4"U|) . 

To describe a rotational band, one should keep the corresponding configuration fixed for 
a set of increasing values of rotational frequency Q. Because of the considerable i^-mixing in 
TAC solutions and the high level density in the 3-dimensional calculations, it may occur that 
the configuration is changing with the iteration and also the different rotational frequency. 
To constrain the specific configuration we are interested in, one can adopt the following 
prescription: starting from the Dirac level \ipi(Q n )) blocked for Q = Q n , one can block for 
Q = fi n +i the level \i/jj(fl n+ i)) which maximizes the overlap (ijji(Q n )\ijjj(Q n+ i)) , i.e., 



(48) 



For infinitesimal step sizes, this condition corresponds to the so-called parallel transport 121 
122 1 . 



III. CELEBRATED MAGNETIC ROTATION IN 198 PB 



A. Brief historical overview 



The observation of the A = 1 rotational-like structures in neutron deficient Pb nuclei 



in the early 1990s opened a new era for magnetic rotation |8HlO|. Lon g cascades of Ml 
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123 



124]. With 



transitions were firstly observed in Pb nuclei in the early 1990's 
improved detector techniques and lots of efforts, in 1997 the lifetime measurements based 
on the Doppler-shift attenuation method (DSAM) for four Ml-bands in the nuclei 198 Pb and 
199 Pb provided a clear evidence for magnetic rotation 29J]. Subsequently, another experiment 



using the recoil distance method (RDM) in 198 Pb together w ith t he results of the DSAM 



125] . Many works along this 



experiment provided further support to the shears mechanism 
line for almost 10 years have been devoted to the magnetic rotation in Pb isotopes which 
becomes a classic example and has received wide attention from theoretical and experimental 
aspects. Naturally, it is a good test ground for various theory aiming at the description of 
the MR. 

In the framework of the pairing plus quadrupole model, the tilted axis cranking ap- 
proach , |6l| has been used to describe the magnetic rotation from the very begin- 



ning and reproduce the experimental reduced B (Ml) values of magnetic dipole bands in 
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198,1991 



model 



3 ve ry well 29[. In addition, the shell model 126] and the many-particles-plus-rotor 



1271 ] have also been applied to investigate magnetic rotation in the Pb region. 



With its many success in describing nuclear phenomena in stable as well as in exotic 



nuclei 



69 



7 If , the CDFT has been generalized to the cranking CDFT |89| and the tilted axis 



cranking CDFT 92|, |93|, |98[ . In particular, the newly developed tilted axis cranking covariant 



98 



includes 



density functional theory (TAC-CDFT) based on point-coupling interactions 
significant improvements by implanting the simplex symmetry and reduces computation 
time, which makes it possible to investigate the MR in the heaviest mass region observed so 
far. 

As the magnetic r otat ion in Pb isotopes is typical and famous examples, in the following 
the band 1 in 198 Pb 123| will be chosen as an example to demonstrate the relativistic self- 
consistent description of the MR characteristics. The calculations J99| have been carried 
out with the covariant point-coupling density functional PC-PK1 [86J. The detailed and 
numerical techniques can be found in Ref. 99[ and references therein. 



B. Single particle Routhian and configuration 



In Ref. 99], TAG based CDFT calculation ha s be en performed with the proton con- 
figuration 7r[s^, 2 2 /i9/2ii3/2]ll~ suggested as in Ref. 123j . The single particle Routhians for 
the neutrons in the nucleus 198 Pb are respectively shown as a function of the rotational 
frequency Q for the two configurations AE11 and ABCE11 in Fig. [51 The positive parity 
levels belonging to the 1^13/2 orbit are given by full black curves and the neutron levels with 
negative parity (pf) are indicated by dashed red curves. In band 1, a backbending phe- 
nomenon has been observed caused by the alignment of a pair of 213/2 neutrons. Before the 
backbending, the neutron configurations ^13/0 ( fa) - . 1 ] h. as b. een assigned to ba nd 1 in 198 Pb. 



After the backbending, it becomes as i / [ii^/ 2 (fp) *] 31 



128) . As in Ref. [124j], a short hand 



13/2* 

notation is used for these configurations: A, B, C and D denote z/z 13 / 2 holes with positive 
parity and by E denotes a neutron hole with negative parity (mainly of f§/2 and P3/2 origin). 
The proton configuration n[s^? 2 hg/2ii3/2}^~ is abbreviated by its spin number 11. Therefore 
the neutron configurations ^[^^(/p) -1 ] is referred as AE11 before the backbending, and 
v[ii S / 2 (/p) -1 ] as ABCE11 after the backbending. 
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FIG. 5: (Color online) Single particle routhians for the neutrons in 198 Pb as a function of the 
rotational frequency based on the configurations AE11 and ABCE11. The blue dots indicate the 
occupied levels at O = and the green dots indicate the occupied levels at the band heads with 



the configuration AE11 
text. Taken from Ref. 



eft panel) and ABCE11 (right panel). Further details are given in the 



C. Energy spectra 



Th e ca lculated energy spectra of the band 1 in Pb are show in comparison with the 



data 



1281 ] in Fig. [6j For certain regions of angular momenta, the calculated values are 



missing, as for instance, / = 19 — 21ft in band 1 in 198 Pb. As discussed in Ref. 93[, it is due to 
the level crossing connected with the backbending phenomenon, and no converged solutions 
could be found for these angular momentum values. It can be seen that the TAC-CDFT 
calculations reproduce well the experimental energies but underestimate the particle-hole 
excitation energies at the band hea d of t he configurations ABCE11 in 198 Pb. In comparison 
with the PQTAC calculations fll . 



129J, these underestimations can be explained by the 
pairing correlations and will be further investigated in the future. At the moment, these 
underestimations are compensated by choosing different references for the configurations 
involved. 

The experimental rotational frequency f2 exp is extracted from the energy spectra by the 
relation 

ft^cxp « — = -[£ 7 (I + 1^I) + E 7 (I^I- 1)]. (49) 
In Fig. [7J the calculated total angular momenta of the band 1 in 198 Pb as functions of the 
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FIG. 6: (Color online) Energy spectra in the TAC-CDFT calculations compared with the data 128 1 
for band 1 in 198 Pb. The energies at I = I5h is taken as references for the band 1 in 198 Pb. Energies 
for the configurations ABCE11* in 198 Pb are renormalized to the energies at I = 22H. Taken from 
Ref. y. 



rotational frequ ency are shown in comparison with the experimental data 128| and the 
PQTAC results [l29j. It is found that both the TAC-CDFT and the PQTAC results agree 
well with the experimental data. This shows that the TAC calculations can reproduce the 
relative changes of the moment of inertia within the different bands rather well. The TAC 
calculations support that the backbendings arise through an excitation of a neutron-hole 
pair in the i\s/2 shell, i.e. by the transitions in the configurations AE11— t-ABCEII in band 
1 of 198 Pb. In detail, before the backbending the spins values found in the TAC-CDFT 
and PQTAC models differ from experimental values up to 2h. After the backbending, the 
PQTAC result for the band 1 in 198 Pb is nearly 3h larger than the experimental values and 
the TAC-CDFT results. Comparing with the experimental values in Fig. [7J the appearance 
of backbending is clearly seen for each band. 



D. Deformation evolution 



The advantage of the TAC-CDFT calculations is that the nuclear shape and deformation 
can be obtained self-consistently and automatically as outputs. The deformation param- 
eters and 7 as well as their evolutions for band 1 in 198 Pb obtained in the TAC-CDFT 
calculations are compared with the PQTAC results 129| in Fig. |HJ In the TAC-CDFT calcu- 
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FIG. 7: (Color online) Angular mome nta a s functions of the rotation al fre quency in the TAC-CDFT 



calculations compared with the data 



1281 ] and the PQTAC results 



129[ | for band 1 in 198 Pb. The 



configurations with "(PQ)" denote the corresponding results of PQTAC calculations. Taken from 



Ref. 



99(]. 
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FIG. 8: (Color online) Deformation parameters /3 (left panels) and 7 (right panels) as functio ns of 



the rotational frequency in the TAC-CDFT calculations compared with the PQTAC results 
for band 1 in 198 Pb. Taken from Ref. 



129] 



lations, the quadrupole deformations are around /3 = 0.15 and remain almost constant. The 
PQTAC calculations produce the same tendency with slightly smaller deformations around 
(3 = 0.11. Meanwhile, the deformation 7 vary between 47° and 59° which means small 
triaxiality close to oblate axial sy mme try in the TAC-CDFT calculations. This is consistent 
with the PQTAC results of Ref. [129 ]. 
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FIG. 9: (Color online) Composition of the total angular momentum at both the minimum and 
the maximum rotational frequencies in TAC-CDFT calculations for band 1 in 198 Pb. Left (right) 



panel is result for the rotation before (after) backbending. Taken from Ref. 
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E. Shears mechanism 

One of the fundamental and important characteristic in magnetic rotation is the shears 
mechanism. In Fig. |9l the proton and neutron angular momentum vectors J n and J v as 
well as the total angular momentum vectors J tot = J n + J u at both the minimum and the 
maximum rotational frequencies in TAC-CDFT calculations for the band 1 in 198 Pb are 
shown. The proton and neutron angular momenta J n and J v are defined as 



Z 

E 

P =i 



(p\J\p), 



(Ju) 



N 
n=l 



n\j\n), 



(50) 



where the sum runs over all the proton (or neutron) levels occupied in the cranking wave 
function in the intrinsic system. 

For the magnetic dipole bands in 198 Pb, the contributions to the angular momenta come 
mainly from the high j orbitals, i.e., the 213/2 neutron (s) as well as /19/2 and 213/2 protons. 
At the band head, the proton particles excited across the closed Z = 82 shell gap into the 
hg/2 an d H3/2 orbitals contribute to the proton angular momentum along the short axis, and 
the neutron hole(s) at the upper end of the 213/2 shell contribute to the neutron angular 
momentum along the long axis. By comparing the upper panels (before backbending) with 
the lower ones (after backbending) in Fig. [91 one finds that after the backbending the neutron 
angular momentum vectors are considerably larger, because they contain the contributions of 
an aligned pair of i 13 / 2 neutron holes. Therefore, the proton and neutron angular momentum 
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FIG. 10: (Color online) B{M\) values as functions of the total angul ar mo mentum in the TAG 



CDFT calculations compared with the data and the PQTAC results 

□ 

Circles an d sq uares denote experimental data from DSAM |29j and RDM 



from Ref. 



99] 



129H for band 1 in 198 Pb. 



1251 ]. respectively. Taken 



vectors form the two blades of the shears. With the frequency increasing, the two blades 
move toward each other and the closing of shears increase the angular momentum, while the 
direction of the total angular momentum stays nearly unchanged. In such a way, the shears 
mechanism is clearly observed. 



F. Electric and magnetic transitions 



A typical characteristic of magnetic rotation is the strongly enhanced Ml transition 
probabilities which decrease with the spin. In Fig. [TOj the calculated B(M1) values as 



functions of the tota 
with the data 129 



angular momentum for the b and 1 in 198 Pb are shown in comparison 
125] and the PQTAC results [129|. The TAC-CDFT calculations reproduce 
the decrease of the observed B(M1) values with increasing spin. However, as observed 
already in earlier investigation 92j, |98j, the absolute values show discrepancies. As shown in 
Fig. [TOl one has to attenuate the results by a factor 0.3 in order to reproduce the absolute 



B(M1) values. The same factor has been used in Refs. 92j, |98j|. So far the origin of this 



attenuation factor is not understood in detail. As discussed in Ref. 



there are however 



several reasons: (a) Pairing correlations strongly affect the deformation and the levels in the 
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FIG. 11: (Color online) B{E2) values as functions of the total angular momentum in t he T AG 



CDFT calculations compared with the data from DSAM 



band 1. Taken from Ref. 



99] 



29(] and the PQTAC results 
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for 



neighborhood of the Fermi surface. This causes a strong reduction for the 5 (Ml) values 
with major contributions from the valence particles or holes, (b) The c oupling to complex 
configurations such as particle vibrational coupling (Arima-Horie effect 
all cases to a quenching of the B(M1) values for neutron configurations 



; 11301. Ill 

mm 



131] ) leads in 



133] . (c) Meson 

exchange curr ents and higher corrections also cause a reduction of the effective ^-factors for 



the neutrons 134j - |l36j ]. However, it is not the absolute B(M1) values, which characterize 
the shear bands, but rather the behavior of these values with increasing angular momentum. 
On the other side, the absolute values of PQTAC results agree with the observed 5 (Ml) 
data and the attenuated TAC-CDFT results. However, they show a sharper decreasing trend 
as compared with the TAC-CDFT calculations. The agreement between the calculated and 
experimental B (Ml) values and their trend shows a convincing confirmation of the shears 
mechanism. 

In contrast to the enhanced Ml transitions, the E2 transitions are weak for magnetic 
rotational bands. In Fig. dU the calculated B(E2) values as functions of the total angular 
momentum ar e sho wn and compared with the DSAM-data of Ref. 29[ and the PQTAC 



results of Ref. (l29| for band 1 in 198 Pb. The B(E2) values in the TAC-CDFT calculations 
are in reasonable agreement with the data and show a roughly constant trend. This is 
consistent with the nearly constant quadrupole deformation in each configuration calculated. 
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Compared to the PQTAC results, the TAC-CDFT calculations predict larger B(E2) values, 
in accordance with the larger deformations shown in Fig. [7J 



IV. MAGNETIC ROTATION IN OTHER MASS REGION 



A. A ~ 60 mass region 



To date, the magnetic dipole bands observed have been summarized in the nuclear chart 
in Fig. El The recent observations in 58 Fe 44] and 60 Ni {45 1 is identified as the lightest mass 
region to exhibit magnetic rotation phenomenon and have extended the observed MR mass 
region to A ~ 60 mass region. 

In Ref. 98j, TAC based CDFT calculation with PC-PK1 86| has been performed for four 



45]. Same as 



magnetic dipole bands, denoted as M-l, M-2, M-3, and M-4, reported in 60 Ni 
in Ref. 45|, the bands M-l and M-4 are suggested to be built from the same type of config- 
urations, i.e., n[(lf 7 / 2 y l {fp) 1 } ® v[{lg 9 /2) l {fpf]- For the bands M-2 and M-3, TAC based 
CDFT calculation clearly indicates that they are respectively built from the configuration 
7r[(l/7/2)" 1 (^9/2) 1 ]®^[(1^9/ 2 ) 1 (/p) 3 ] ^&ii[(lf\ /2 y\f P y}®u[(lg, /2 )\f P ) 2 }. For simplicity, 
the above configurations are referred as Configl, Config2, and Config3, respectively. 
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FIG. 12: (Color online) Energy spectra obtained from the TAC-CDFT calculations in comparison 
with the available data for bands M-l and M-4 (left panel), M-2 (middle panel), as well as M-3 
(right panel). The energies at / = 9h, I = llh, and / = 15h are taken as references in the left, 
middle and right panels, respectively. The evolutions of the nuclear shape (0,^) for bands M-l, 
M-2, and M-3 are also illustrated with the schematic pictures. Taken from Ref. 



The calculated energy spectra are shown in compared with the available data for the 
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bands M-l and M-4 (left panel), M-2 (middle panel), as well as M-3 (right panel) for 60 Ni in 
Fig. [T2j In general, the experimental energies of the bands M-l, M-2, and M-3 are reproduced 
very well by these TAC-CDFT calculations. However, the assigned configuration for each 
of these bands could not be followed in the calculations up to the highest spin observed, 
i.e., convergent results could be obtained only up to ~ 14ft for Configl, ~ 15ft for Config2, 



and ~ 16ft 
evolution 



'or Config3. These are connected with the configuration change and shape 



With the increase of the rotational frequency, the configurations, 7r[(l/ 7 / 2 ) 1 (fp) 1 ] <8> 
v[(lg B/2 ) 1 (fp)*(lf 7/2 )- 1 ] (Configl*) and7r[(l/ 7/2 )- 2 (/p) 2 ]®z/[(l^ 9/2 ) 2 (/p) 3 (l/ 7/2 )- 1 ] (Con- 
fig3*) will strongly compete with Configl and Config3, respectively. In other words, one 
observes a neutron pair broken in the / 7 / 2 shell at / = 15ft in band M-l, and the excitation 
of a unpaired proton from the / 7 / 2 shell to the fp orbital and a neutron pair broken in the 
/ 7 / 2 shell at / = 16ft in the band M-3. 

The shape evolutions of bands M-l, M-2, and M-3 are also shown in Fig. [T2] It is 
interesting to note that the nucleus changes its shape from prolate-like to oblate-like with the 
frequency in Configl, Config2, and Config3, and comes back to a prolate-like deformation 
with the configuration changing from Configl to Configl*, and Config3 to Config3*. In 
particular, the nucleus with Config3* has a relatively large deformation (/3 ~ 0.4) with axial 
symmetry. 

For the magnetic rotation bands in 60 Ni, the proton and neutron angular momentum 
vectors and J u defined in Eq. ( 1501) as well as the total angular momentum vector «7 tot at 
both the bandhead and the maximum rotational frequency in the TAC-CDFT calculations 
with the configurations of Configl, Config2, Config3, and Config3*, are shown in Fig. [13j 

For the bands built on Configl, Config2, and Config3, the contributions to the angular 
momenta come mainly from the high j orbitals, i.e., the gg/ 2 neutron(s) and the / 7 / 2 proton. 
At the bandhead, the neutron particle(s) filling the bottom of the gg/ 2 shell mainly contribute 
to the neutron angular momentum along the x-axis, and the proton hole at the upper end of 
the fj/2 shell mainly contributes to the proton angular momentum along the z-axis. They 
form the two blades of the shears. As the frequency increases, the two blades move toward 
each other to provide larger angular momentum, while the direction of the total angular 
momentum stays nearly unchanged. In this way, the shears mechanism is clearly seen. 

One should notice that the proton particle in the g$/2 orbital also give substantial con- 
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FIG. 13: (Color online) Composition of the total angular momentum at both the bandhead and the 
maximum rotational frequency in the TAC-CDFT calculations with the configurations of Configl, 
Config2, Config3, and Config3*. Taken from Ref. [98]. 



tributions to the proton angular momentum in the case of Config2. As a result, J n has not 
only a large J z component but also a substantial J x component even at the bandhead. Thus, 
the shears angle G, the angle between J n and J u , is not as large as those of Configl and 
Config3, and decreases only by a small amount with increasing rotational frequency. 

For the Config3*, as the two proton holes in the f?/2 orbital are paired, the proton angular 
momentum comes mainly from the particles in the fp shell, which aligns along the x-axis. 
The neutron hole in the fan orbital gives substantial contributions to the neutron angular 
momentum, which leads to a large J z component. Higher spin states in the band are created 
by aligning the neutron angular momentum towards the x-axis. Considering the large axially 
symmetric prolate deformation as shown in Fig. [T2J the mechanism of producing higher 
spin states with Config3* is electric rotation rather than magnetic rotation. Therefore, a 
transition from magnetic rotation to electric rotation is observed in Config3*. 

Recently, the high spin structure in 58 Fe has been investigated by heavy-ion induced 



fusion-evaporation reactions at Gammasphere |44j. The magnetic rotational bands observed 
have been interpreted with the TAC-CDFT, which concludes that 58 Fe is the lightest nucleus 
exhibiting magnetic rotation. 
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B. A ~ 140 mass region 



For A ~ 140 mass region, nhn/2 particles combined with vh\xj2 holes satisfy the high-j 
configurations for magnetic rotation. In fact, the magnetic rotation in A ~ 140 mass region 
have been identified in Tb, Gd, Eu, Sm, Pm, Nd, Pr, Ce, La, Ba, Cs, Xe, Te and Dy 
isotopes (see 



311 ] and the references therein). 
For 142 Gd, five bands denoted as DB1, DB2, DB3, DB4 and DB5, have been observed 
and four of them have been interpreted as magnetic rotation bands with the con figuration s 



ixh\ l/2 



^11/2' 7r ^n/2 



z//i n 4 /2 , nh\ 1/2 <g> 7rg 7 ^ 2 uh^ /2 , and 7rh\ 1/2 <g> Kg 7 f 2 vh 1 { /2 



'huh-, 



1371 . 



138|, 



respectively. In Ref. 66[, the shears mechanism and the spectrum of the band DB1 in 
i been investigated by the tilted axis cranking Skyrme Hartree-Fock method. In 
TAC-CDFT with the meson-exchange interaction PK1 



the magnetic rotational band DB1 in Gd based on the configuration tt/i^ 2 



1391 ] has been applied for 



v\l j2 . For 



simplicity, the configurations ith^^ (g> vh xl , 2 and ir[hl 1 ^ 2 g 7 ^ 2 d 
Configl and Configl*, respectively. 
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5/2J 
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FIG. 14: (Color online) The energy as a function of the total spin in TAC-CDFT calculation [93] 
for Configl (full red dots) and Configl * (gr een triangles up) in comparison with the data for the 



band DB1 in 142 Gd (filled back circles) 



1381 ]. The non-relativistic SHF result (blue triangles down 



661 ] is also included. 



The energies as functions of the total angular momentum for the Configl calculation 
(full red dots) and the Configl* calculation (green triangles ups) for the magnetic dipole 
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band D B1 in 142 Gd are shown in Fig. [TH They are compared with the available data (filled 
circles) 1381 ] and the non-relativistic SHF results from Ref. 66j (triangles down). As no li nk 



to the ground state is observed, the recommended band head spin / = 16 ft of Ref. 138j is 
adopted, and the energy at I = 16 h are taken as a reference for both the RMF values and 
the nonrelativistic SHF calculation. In general, the energies in RMF calculations achieve 
better agreement with the data as compared with the SHF results. 
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FIG. 15: (Color online) The total angular momentum as a function of the rotational frequency 



in TAC-CDFT calculation 



for the configurations vr/i^ , 2 (g) vh-^^ (Configl, open circles) and 



/2 



^ [^ll/2#7/2^5/3. 

(filled circles) 



1381 ] 



^11/2 (Configl*, triangles up) in comparison with the data for DB1 in Gd 



The total angular momenta are shown as functions of the rotational frequency and re- 
sults of Configl (full red dots) an d Co nfigl* (green triangles up) calculations are compared 
with data (filled circles) of Ref. 138j in Fig. [15j The total angular momenta in Configl 
linearly increase with the rotational frequency and agree with the data till HQ = 0.40 MeV. 
After Ml > 0.40 MeV, a up-bending is observed for the data. This up-bending can not be 
reproduced by the smooth behavior in either Configl or Configl* calculations. 



C. A ~ 80 and A ~ 110 mass regions 



For A ~ 80 mass region, 71779/2 particles combined with ugg/2 holes satisfy the high-j 
configurations for magnetic rotation. In fact, the magnetic rotation in A ~ 80 mass region 
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have been identified in Rb, Kr and Br isotopes (see 3l| and the references therein). 

As the first relativistic investigation of the magnetic rotation, three-dimensional cranking 
CDFT has been developed and applied for 84 Rb 92[. The proton configuration is fixed to 
be 7r(p/) 7 (lg 9 /2) 2 with respect to the Z = 28 magic number and ^(lg^) -3 with respect 
to the N = 50 magic number is adopted for the neutron configuration. The signals of the 
shears mechanism, such as the nearly constant tilt angle and the smooth decreases of the 
shears angle and of the B(M1)/B(E2) ratio, are well reproduced. The detailed discussion 



can be seen in Ref. 



92] . Because of the numerical complexity, so far, the three-dimensional 



92). 



cranking CDFT has been applied only for the magnetic rotation in Rb 

For A ~ 110 mass region, irhn/2 particles combined with ugg/2 holes satisfy the high-j 
configurations for magnetic rotation. In fact, the magnetic rotation in A ~ 110 mass region 
have been identified in Cd, In, Sn, Sb, and Te isotopes (see [3l| and the references therein). 
Using the recently developed TAC-CDFT based on point-coupling interactions, the magneti c 
rotation bands in 113 > 114 In in A ~ 110 mass region are well reproduced successfully 140lll41 |. 

In order to explore the MR in the lightest mass region, the TAC-C DFT calculation has 
been performed for 22 F with the configuration <S> v d~^, 2 in Ref. 1421 ] . The possible 

existence of magnetic rotation is suggested for 22 F via investigating the spectra, the relation 
between the rotational frequency and the angular momentum, the electromagnetic transi- 
tion probabilities B(M1) and B(E2) together with the shears mechanism characteristic of 
magnetic rotation. 



V. ANTIMAGNETIC ROTATION 



A. Concept 

As mentioned in the introduction, in analogy with an antiferromagnet, antimagnetic 
rotation" (AMR) [9] is predicted to occur in some specific nearly spherical nuclei, in which the 
subsystems of valence protons (neutrons) are aligned back to back in opposite directions and 
nearly perpendicular to the orientation of the total spin of the valence neutrons (protons). 
Such arrangement of the proton and neutron angular momenta also breaks the rotational 
symmetry in these nearly spherical nuclei and causes excitations with rotational character 
on top of this bandhead as antimagnetic rotation 
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106/ 



To date, antimagnetic rotations have been reported in Cd isotopes inclu ding: 105 Cd 48], 
Cd |42|, 108 Cd 49, 50), and 110 Cd 52]. The other candidates include 109 Cd 51], 100 Pd 54], 



and 144 Dy [55(. In order to apply the CDFT for antimagnetic rotation phenomenon, the 
newly developed TAC-CDFT based on point-coupling interactions are used to investigate 
anti magn etic rotation (AMR) in 105 Cd in a fully self-consistent and microscopic way in 



Rcf. 



100 |. 



B. Energy spectrum 




FIG. 16: (color online) Energy (upper panel) and rotational frequency (lower panel) as functions 
of the total angular momentum. The fully self-consistent solution (solid lines) and that neglecting 



polarization (dashed lines) are compared with the data 
is taken as reference in the upper panel. Taken from Ref 



4£l (solid dots). The energy at / = 23/2ft 



10C|. 



For the TAC-CDFT calculation 



1001 ] based on point-coupling interactions PC-PK1 86] 



for AMR band in Cd, the odd neutron occupies the lowest level in the hn/2 shell and 
the remaining nucleons are treated self-consistently by filling the orbitals according to their 
energy from the bottom of the well. This automatically leads to the configuration for AMR. 



In 

lines) 



ig 



[TBI the TAC-CDFT calculated energy and the rotational frequency (solid 



XOOf ] are compared with data 481 ] ■ In the upper panel it can be clearly seen that, 



apart from the bandhead, the experimental energies are reproduced excellently by the present 
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self-consistent calculations. In the lower panel it is found that the calculated total angular 
momenta agree well with the data and increase almost linearly with increasing frequency. 
This indicates that the moment of inertia is nearly constant and well reproduced by the 
present calculations. 



C. Two shearslike mechanism 
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FIG. 17: (color online) Angular momentum vectors of neutrons J v and the two gg/ 2 proton holes 
at bot h the bandhead (Ml = 0.3 MeV) and the maximum rotational frequency. Taken from 



Ref. 



10C]. 



In order to demonstrate the two shearslike mechanism in 105 Cd in Fig.[T71 there was shown 
that both at the bandhead and at the maximum rotational frequency the angular momentum 
vectors of the two g$ /2 proton-holes j n and of the neutrons J u = 3^ where n runs over 
all the occupied neutron levels. At the bandhead, the two proton angular momentum vectors 
j n are pointing opposite to each other and are nearly perpendicular to the vector J u . They 
form the blades of the two shears. With increasing Q the gradual alignment of the vectors 
jV of the two gg/2 proton holes toward the vector J v generates angular momentum while 
the direction of the total angular momentum stays unchanged. This leads to the closing of 
the two shears. The two shearslike mechanism can thus be clearly seen, and it is consistent 
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with the previous works 
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FIG. 18: (color online) Alignment of the valence neutrons (left side) and proton holes (right si de) a t 



both the bandhead (Ml = 0.3 MeV) and the maximum rotational frequency. Taken from Ref. 



10C]. 



In a microscopic calculation, there is no inert core and all the energy and angular momen- 
tum comes from the particles. The contributions of the valence neutrons and proton holes 
to the angular momentum J x at both the bandhead and the maximum rotational frequency 
are shown in Fig. [18j It is found that the contributions come mainly from high-j orbitals, 
i.e., from gg/2 proton holes as well as from hu/2 and g^/2 neutrons. In order to provide a 
simple picture which can be compared with the core angular momentum given in Ref. {9J, 
one can estimate the "corelike" angular momentum in the present framework by excluding 
the contributions of three valence neutrons, shown in the left three columns in Fig. 2, from 
the total neutron angular momentum. It is found that the "core" contributes about 3 h 
when the frequency Q increases from the bandhead to the maximum value. 

For the protons, only the two holes in the gg/2 shell contribute. As shown in Fig. [TTl 
they cancel each other in the z direction giving non-negligible contributions to the angular 
momentum along x axis even at the bandhead. With growing frequency, the proton angular 
momentum in x direction increases because of the alignment of the two proton hole blades. 
For the neutrons, on the other hand, we have only contributions above the A^ = 50 shell. 
One neutron sits in the hu/2 orbit and the other six are, because of considerable mixing, 
distributed over the gj/2 and d^/2 orbitals. As Q grows, the contributions of the aligned 
orbitals with j x = 11/2 and 9/2 do not change much and the increase in angular momentum 
is generated mostly by the alignment of orbitals with low j x values. This microscopic calcu- 
lation shows that the interpretation given in Ref. 4{| is only partially justified: we clearly 
have two proton holes in the g 9 / 2 and one neutron particle in the h u / 2 orbit, but, due to the 
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mixing of orbits with lower j values the other neutrons are distributed over several subshells 
above the N = 50 core and the increasing angular momentum results from the alignment of 
the proton holes and the mixing within the neutron orbitals. Because of this str ong mixing 



between the neutrons, a core needed for the phenomenological model in Ref. 
really be defined. 



10] cannot 



D. Electric transition probability and deformation 
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FIG. 19: (color online) B{E2) values as a function of the angular momentum. Solutions with 



(solid line) and without (dashed line) polarization are compared with data 



48[ (solid dots). Inset: 



Deformations (3 and 7 drive n by the increasing rotational frequency whose direction is indicated 



by arrows. Taken from Ref. 



10C] 



AMR is characterized by weak E2 transitions decreasing with increasing spin. In Fig. [T9~| 
the calculated B(E2) values (solid lines) are compared with the available data 48|. It 
is found that the resulting B{E2) values are very small (< 0.14 e 2 b 2 ) and in very good 
agreement with the data. Furthermore, the fact that the B(E2) values decrease with spin 
is in agreement with the interpretation by the two shearslike mechanism. 

The decrease of the B(E2) values can be understood by the changes in the nuclear defor- 
mation. As shown in the inset of Fig. [19], with increasing frequency, the nucleus undergoes 
a rapid decrease of deformation from 0.2 to 0.14 at a small and near-constant triaxiality 
(7 < 9°). As usual, it is found that the deformation of the charge distribution, responsible 
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for the B(E2) values, changes in a similar manner. Therefore, one can conclude that the 
alignment of the proton and neutron angular momenta, i.e., the two shearslike mechanism, 
is accompanied by a transition from prolate towards nearly spherical shape. 



E. Core polarization 

In order to investigate the importance of the polarization effects induced by the two 



proton holes, which is taken into account fully in the TAC-C DF 



calculation without polarization has been carried out in R ef 
frequency the 107 Sn core is firstly calculated in Ref. 



calculation, an additional 



100J. For this purpose, at each 



100], where the two proton holes are 



filled. This results in a filled nearly spherical g$/2 shell. In the next step, the self-consistency 
(dashed lines in Figs. [TBI and [T9"]) is neglected and the band in 105 Cd is calculated, in which two 
protons are removed from the gg/2 shell using, at each frequency, the corresponding nearly 
spherical potentials S and obtained in the calculations of the 107 Sn core. As shown in 
the upper panel of Fig. [IBJ the energy is reduced only slightly by polarization in the lower 
part of the spectrum. At the same time the relation between angular velocity and angular 
momentum is considerably changed in the lower part of Fig. [161 Without polarization, 
a much smaller frequency Q is necessary to reach the same angular momentum as with 
polarization. In addition, without polarization, there is a maximal angular momentum of 
roughly 17 H. Higher values cannot be reached even at rather high frequencies. 

This behavior can be well understood from the evolution of the deformation parameters 
shown in the inset of Fig. [19j Without polarization we use at each frequency the potentials 
of the core-nucleus 107 Sn, where the deformation is relatively small. Angular momentum 
can only be produced by alignment of neutron particles along the rotational axis leading at 
Hfl = 0.5 MeV to an oblate shape with a rotation around the symmetry axis. Removing two 
protons would lead, if polarization is taken into account, to a larger prolate configuration 
with lower energy and hindering alignment. Therefore, in the lower part of the band, it is 
easier to produce angular momentum without polarization, where the deformation is small. 
On the other side the oblate deformation keeps the high-j proton holes in the j x = ±1/2 
orbitals of the gg/2 shell pairwise occupied and hinders their alignment. Above hQ = 0.5 MeV 
we reach the maximum angular momentum of the neutron configuration. With polarization, 
because of the prolate deformation, we have much more mixing and therefore can reach 
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larger angular momentum by aligning the protons. This can also be seen in Fig. [19] that 
the E2 transitions do survive with polarization when Ml > 0.5 MeV. Therefore, it is of 
importance to emphasize that polarization effects play a very important role in the self- 
consistent microscopic description of AMR bands. 

VI. SUMMARY AND PERSPECTIVES 

In the past decades, the rotational-like sequences in near-spherical or weakly deformed 
nuclei have attracted significant attentions. This phenomenon, known as magnetic rotation, 
has been extensively explored experimentally and theoretically. With its many success in 
describing nuclear phenomena in stable as well as in exotic nuclei, the CDFT has been 
generalized to the cranking CDFT, the tilted axis cranking CDFT and three-dimensional 
cranking CDFT and applied for electric and magnetic rotations all over the nuclear chart. In 
particular, the newly developed TAC-CDFT based on point-coupling interactions includes 
significant improvements and reduces computation time, which makes it possible to perform 
systematic investigation. 

This review provides an overview of the experiential and theoretical status of MR and 
AMR, a sketch of the CDFT as well as TAC-CDFT based on point-coupling interactions, 
and followed by the summary of the TAC-CDFT descriptions for the MR and AMR bands. 

The shears bands in the nuclei 198 Pb is the most well-known example of MR. The TAC- 
CDFT calculated energy spectra, the relation between spin and rotational frequency, the 
deformation parameters and reduced Ml and E2 transition probabilities are discussed and 
compared with data. By choosing 60 Ni, 84 Rb, and 142 Gd as examples, the success of the 
TAC-CDFT in describing MR for A ~ 60, A ~ 80, and A ~ 140 mass regions has been 
demonstrated. By reproducing the experimental energy spectrum and the B(E2) values 
combining with the examination of the core polarization and deformation evolution, an 
antimagnetic rotation band in 105 Cd has been investigated in a fully self-consistent and 
microscopic way by the TAC-CDFT. The two shearslike mechanism in AMR has been clearly 
illustrated. 

One should note that the pairing correlations have been neglected in the calculations 
presented in this review. In fact, the pairing correlations are negligible for many MR and 
AMR bands in which the high j oribitals near the Fermi surface are blocked and thus reduce 
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the pairing effects. However, in some specific MR and AMR bands, the low j orbitals 
also play important role in the rotation excitations. Therefore, it is very important to 
investigate the pairing effects on the MR and AMR in the framework of the tilted axis 
cranking covariant density functional theory. Due to the time-reversal symmetry broken, 
the simple BCS method is not valid and one should resort to the Bogoliubov method. In 
such case, the treatment of particle projection is usually necessary since the particle number 
conservation is violated in the Bogoliubov transformation. 



There is another way known as Particle Number Conservation (PNC) method 143l Il44j . 



which can treat the pairi ng correl ation exactly and has been implemented successfully in the 



cranking Nilsson model [144H146I ]. Similar method, the Shell- model-like Approach (S LAP) , 



1471 . 



has also been applied to treat the pairing correlations in the framework of CDFT 
The SLAP has the advantage that the paring correlation is treated exactly and thus the 
particle number conservation is not violated in the calculations. Therefore, it would be very 
interesting to introduce this method to the present tilted axis cranking covariant density 
functional theory to investigate the pairing effects in the tilted axis cranking calculations. 

During the past decades, the three-dimensional cranking calculation has attracted wide 
attentions since it is connected with the novel phenomenon of chiral rotation in nuclei. 
Therefore, it is also a very interesting to extend the tilted axis cranking covariant density 
functional theory to three-dimensional cranking case with the point-coupling interaction. 
Similar work has been performed with the meson-exchange interactions. However, efforts 
should be made to simplify its numerical complicity to investigate the chiral rotation. 
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